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Abstract 



A local, second-order (truncated) approximation is applied to the Hub- 
bard model in three dimensions. Lowering the temperature, at half-filling, 
the paramagnetic ground state becomes unstable towards the formation of a 
commensurate spin-density-wave (SDW) state (antiferromagnetism) and suf- 
ficiently far away from half-filling towards the formation of incommensurate 
SDW states. The incommensurate-ordering wavevector does not deviate much 
from the commensurate one, which is in accord with the experimental data 
for the SDW in chromium alloys. 
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I. INTRODUCTION 



Pure Cr has a spin-density-wave (SDW) ground state with an ordering wavevector which 
is incommensurate with the underlying lattice§. Adding electrons to Cr (by alloying with 
Mn) makes the magnetic order commensurate with the lattice and increases the transition 
temperature, whereas removing electrons from the system (by alloying with V) drives the 
magnetic order to a more incommensurate one, and decreases the transition temperature, 
eventually to zero (Fig. 5). Penni found qualitatively the same behavior for the ground 
state of a single-band Hubbard model within a mean-field approximation. The mean-field 
approximation leads to the usual Stoner criterion for the instability of the paramagnetic state 
to the formation of a SDW state. A first-order perturbation expansion of the self-energy 
in terms of the interaction parameter leads to the same criterion. In this contribution we 
extend this approximation one step further to a second-order (non-self-consistent) expansion 
for the self-energy which includes the lowest-order quantum fluctuations. But, motivated 
by the case for large spatial dimensions, we ignore the nonlocal (site-non-diagonal) elements 
of the self-energyi. In fact, inclusion of the lowest-order quantum fluctuations has been 
shown to have a dramatic effect on transition temperatures and phase diagrams in the 
large-dimensional limit, and has been shown to agree with the quantum Monte Carlo data 
(QMC) over a wide range of interaction strength^Hl. 

The Hubbard mo defl is perhaps the simplest model which can be used to study the many- 
body aspects of correlated electrons on a lattice. The Hubbard Hamiltonian is written in 
the form 

'H = -t J2 i4,aCj,a + c]^^Ci^a) + U ^i,T^i,i (l) 
<i,j>,cr i 

where c|g.(cj^o-) represents the creation (destruction) operator of an electron in a Wannier 
state of spin cr (cr = ±| or t|), on site i, and rij = c-^^Ci^^ is the electron number operator. 
The first term corresponds to the kinetic energy and describes the hopping of electrons 
between nearest-neighbor sites on a lattice via an overlap integral t. This term gives a 
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tight-binding description of the electrons in a periodic potential forming a single energy 
band ecr(k) = —2t Z^f^i cos fc^a for a simple cubic lattice with lattice constant a . The second 
term corresponds to the Coulomb repulsion between electrons. The long-range Coulomb 
interaction is assumed to be screened in the solid so that only the interaction between two 
electrons on the same site is retained, yielding the additional energy of U when the lattice 
site is doubly occupied. 

The model is specified by three parameters: the strength of the electron interaction U 
(measured relative to t) ; the electron density per spin or electron-filling — J2i,a < 
cl f^Ci^a >, where N is the number of sites in the lattice; and the temperature T. 

In Section II we introduce the formalism and the approximation that we use to form 
the phase diagram. In Section III, the details of the numerical calculations are described. 
Section IV, presents the results for the second-order approximation and compares them to 
the first-order approximation, and the QMC results. A semiquantitative comparison is also 
made with the experimental data for Cr. Conclusions follow in Section V. 

II. FORMALISM 

For a given U and Ue there may exist more than one type of spin order for the ground 
state, each being stable at a different temperature. Here we start from the paramagnetic 
state and find a criterion for the instability towards the formation of a SDW state. To do 
this we couple an external magnetic field to the system and look for singularities in the 
response function (magnetic susceptibility) as we change the model parameters. The total 
Hamiltonian of the system in the presence of the magnetic field hi is 

nh = n-Y.hiSt. (2) 

i 

where Ti. is the Hamiltonian in the absence of the external field [found in Eq. (1)] and i is the 
position of the ith electron with spin — c^io-- The spatial variation of the external 
field hi is chosen to probe the particular expected order for the spins. For example, if we 
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want to examine the instability towards antiferromagnetism we choose hi to be of the same 
magnitude everywhere but of the opposite sign on the two sublattices of the bipartite lattice. 

The static response of the system (or the static spin susceptibility) at temperature T is 
defined as follows: 

X.,-2^^U=o, (3) 
= 2Tj2a^\,=oe-'--^-, (4) 

where Q^-" = J rfre^*'^"'^ < T^Cj o-(r)cj^^(0) >, is Green's function at the Matsubara frequency 
Un = {2n + l)7rT. 

Dyson's equation for Green's function (of the Hamiltonian Tih) is, 

{g-''^)-l = {g^^^)T^-T:if + ahA,, (5) 

where Q^"^'" is the noninteracting {U = 0) Green's function and S^-'^ is the matrix element 
of the proper self-energy. Relation (5) and the derivative of the following identity, 

are employed to find, 

Xij = Xij +2 2^ 2^ crg^^ Qn \h=Q, [7) 

k,k',l,l' n',a,a' 'jyk'U ^"'i 

where Xij = TY.nX7j and the bare susceptibility satisfies = —Y^aGij^Qji"- 
For the paramagnetic state, we write Qj'j'^ = G^j''^ and = , so that, 



k,k',Ll' n' 



We need an explicit form for the self-energy to simplify Eq. (8) further. Motivated by 
the work for large spatial dimensional we use a perturbation expansion for the self-energy 
up to second-order in U/t. It was found that the resulting self-energy is almost local in 
three dimensions!, i.e. Sjj k, J^aSij, so we employ the local approximation Sjj = T^ii5ij. 
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This approximation becomes exact in large spatial dimensions!, and in three dimensions, 
the effect of the nonlocal fluctuations on is around 3% in the weak-coupling limitl. 

Fig. 1 shows the diagrammatic expansion of the local self-energy through second order, 
which includes the Hartree term and the second-order bubble. Evaluating the diagrams 
yields. 



m,—cr 
n,n' 



(9) 



for the second-order local self energy. Substituting this approximation in Eq. (8) and Fourier 
transforming to the reciprocal lattice, gives a Dyson's-like equation for the susceptibility 

X(q) = X°(q) + E X°"(q)rr x"'(q), (10) 

n,n' 

where x(q) is the Fourier transform of Xij the first Brillouin zone and 

= U[l-Ux7oA^^n+n')i (11) 

is the irreducible vertex function, xfoc denotes the local particle-particle susceptibility which 
is given by, 

xroc{^^n) = ^ E ^'■(p)^-^+"(p'), (12) 
and x°(q) is the usual bare particle- hole susceptibility, 

x°(q) = ^E^"(p)^"(p + q)- (13) 

We expect x(qc) to diverge at some temperature Tc and electron-filling riec when the 
system is unstable towards the formation of a SDW at the ordering wavevector q^. So, 
near the transition temperature, we can neglect x°(q) compared to x(q) Eq.(lO). At low 
temperatures, the temperature dependence of xfJc is found to be negligible and hence F""'""' 
becomes independent of n and n' . Thus for low transition temperatures and small t^xfJc; ^® 
find from Eq.(lO) and Eq.(ll), 

^ = /(qc, T„ ne.) - xZiT = 0, n^e), (14) 



as the condition for the transition from the paramagnetic phase to an ordered SDW phase. 
This equation is called the modified Stoner criterioni for the magnetic instability of the 
paramagnetic ground state. If the particle-particle susceptibility is ignored on the right- 
hand-side of Eq. (14), the modified Stoner criterion becomes the usual Stoner criterion; this 
term results from including the second-order graph in the self-energy expansion. 

We must apply one more approximation to calculate the two types of susceptibilities that 
appear in modified Stoner criterion . We use bare Green's functions in Eq.(12) and Eq.(13) 
to calculate these susceptibilities. In other words, we are performing a truncated expansion 
for the self-energy that is not self-consistent. This yields 

/(q) = y [ dPp \- . ^ , (15) 

(27r)^V-^ - e(p) + zcj„ - e(p + q) + /i 

and 

(T = 0) = lim V 

Xiocl^ nm ^^^^^ ^ 

The lowest- order effect of the quantum fluctuations is remarkably simple: just reduce the 
momentum- dependent particle-hole susceptibility by the local particle-particle susceptibility 
before applying the Stoner criterion. 



d^p 



1 

iun - e(p) + /i 



(16) 



III. NUMERICS 

For each value of Tc, qc and f/, the root of Eq.(14) yields the critical filling Uec for the 
SDW order. We use the particle-hole symmetry of the model to find the phase diagram 
only for values of Ue < 0.5; the phase diagram is symmetric around Ug = 0.5. For a fixed 
temperature we can find different roots by changing the ordering wavevector qc. Since the 
modified Stoner criterion of Eq. (14) is valid only in the nonmagnetic region of the phase 
diagram, we accept that value of filling for which an instability occurs first, as we approach 
half-filling, i.e., we search for the qc in the first Brillouin zone which makes the filling minimal 
for a fixed temperature. 
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Calculations of the modified Stoner criterion for finite-size lattices show that the desired 
changes only along the edge of the Brillouin zone as the system is doped away from 
half-filling, i.e. qc = (tt, tt, q^) for the reduced zone that contains the z-axis. We use this 
fact to reduce the triple integral in Eq.(15) to an effectively one-dimensional integral in the 
following way: First note that we can write, 

-If,. 1 



(27r)3 J 2{iuJn + /u) + 2tcos{pz) + 2tcos{pz + Qz) 

(17) 



1 1 

+ - 



_iu>n - e(p) +11 iuJn- e(p + q) + ii_ 

This becomes an integral over the 3D density of states ii Qz — n (i.e. at the zone corner). 
Note however that the dependence on p^ and Py is only through —2tcos{px) — 2tcos{py) 
which allows the p^ and Py integrals to be replaced by an integral over the 2D density of 
states, and produces local Green's functions in 2D: 

X'"(ci) = ^jdp ^ 



271 2{iujn + /x) + 2t cos{pz) + 2t cos{pz + Qz) 

■ [G2D{il^n + 2t COs{pz) + A*) + G2D{il^n + 2t COs{pz + Qz) + A*)] • (18) 

Here Q2d{z) is the local 2D Green's function, 

QMz) = I de^. (19) 
J z — e 

G2d{z) is evaluated with a quadrature algorithm which employs a rational function ex- 
pansion for large \z\, and employs a 512 point Gaussian integration when \z\ is small. So 
we can evaluate x°(q) efficiently by numerically performing the remaining integration over 
Pz- Note that xfj^ is already in the form of a one- dimensional integral over the 3D density 
of states. 



IV. RESULTS 



Fig. 2 shows the resulting phase boundary for the Hubbard model using the modified 
Stoner criterion. At half-filling the transition is always commensurate with the lattice, i.e. 
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q = ^(1,1,1). As the system is doped away from half-filling the transition temperature 
decreases and a point is reached where the transition becomes incommensurate with the 
underlying lattice. Finite-size calculations show that the ordering wavevector changes only 
along the Brillouin zone edge as the transition becomes incommensurate, i.e. q = ^(1, 1, 1 — 
6) . The inset in Fig. 2 shows that 6 remains small and eventually stops changing as the 
filling is changed. Note that the ordering wavevector initially changes very rapidly at the 
incommensurate order onset, which is reminiscent of the first-order jump seen in the Cr 
datJ. 

Fig. 3 shows the result of a finite-size calculation for the original Stoner criterion which 
results from a first-order (Hartree) approximation for the self-energy. We see that the 
transition temperature is enhanced by almost a factor of four for each value of U compared 
to the second-order approximation results at half-filling. For a given value of U the change in 
slope at the onset of incommensuration is much smaller than the results of the second-order 
approximation. This change of slope in the second-order approximation is again similar 
to what is seen in the Cr data. Also, as was already shown by Penni, in the first-order 
approximation, for values of U large enough, the ordering wavevector changes first along the 
edge of the zone, but then continues to move toward the zone center as the system is doped 
far away from half-filling. This latter behavior has suggested that the system can become 
ferromagnetic for low electron filling which is ruled out once the quantum fluctuations are 
included. 

Fig. 4 shows a comparison of the transition temperature versus U at half-filling for the 
first and the second-order approximations as well as the result of a quantum Monte Carlo 
(QMC) simulation!. We see that in the case of the second-order approximation, Tc has 
negative curvature for U/t around 10 similar to the QMC result (although the former is 
not a valid approximation for large values of U), whereas in the first-order approximation 
Tc continues to increase with U with a positive curvature. In fact, since the second-order 
approximation is exact as f/ — > 0, it is clear that the QMC data is overestimating by a 
factor of four in the weak-coupling limit. Such criticism has already been raised by comparing 
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to other approximation methodsEj. 

Electronic band structure calculations^ show that the c?-electron concentration for pure 
Cr is 2.28 electrons/atom/spin. In order to map this onto our single-band model, we assume 
that 1/5 of this contributes to our single-band filling, implying rte = 0.456 (there are five 
d-bands in Cr). The Fermi energy of pure Cr is near the edge of a (i-band so that the density 
of states has a large peak there, thus we assume that doping with Mn adds one electron to 
the single band, and doping with V removes one electron from it. This maps the rest of 
the experimental data to the single-band model. We then fit the value of U in the modified 
Stoner criterion, to reproduce the experimental filling for the onset of incommensurate order. 
As shown in Fig. 5, the fit is remarkably good for U = 5.5t. Note that there is only one 
adjustable parameter in this fit. The agreement is remarkable since the lattice structure 
for Cr is actually bcc and because we made such simple assumptions. Furthermore, the 
experimental shape cannot be properly accounted for with a first-order approximation for 
any reasonable value of U. 



V. CONCLUSIONS 

We found a criterion for the instability of the paramagnetic ground state of the 3D 
Hubbard model towards the formation of an ordered SDW state. This criterion is a simple 
modification of the Stoner criterion and is derived by including the second-order graph in 
the expansion of the local self-energy. Hence this new criterion includes the lowest-order 
corrections of the old Stoner criterion due to quantum fluctuations. 

We showed that the magnetic phase boundary between paramagnetic and SDW phases 
of the model changes both quantitatively and qualitatively when we include the effect of 
quantum fluctuations to the usual mean- field (Hartree) result. The quantum fluctuations 
suppress the transition temperatures and exclude the possibility of the formation of ferro- 
magnetism in the 3D Hubbard model. 

We showed that the resulting phase boundary from this modified Stoner criterion is 
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consistent with the experimental data for the SDW in dilute chromium alloys whereas there 
are features in the data that cannot be accounted for by the usual Stoner criterion with 
reasonable values of the interaction strength. We were indeed able to achieve remarkable 
agreement with the experimental data, although this model is too simple to be a realistic one 
for Cr. However, since the modified Stoner criterion is a simple and efficient approximation, 
we hope it can be applied to more realistic models for Cr and other transition metals. 
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FIGURES 

-G 




FIG. 1. The local self-energy Ejj in the second-order approximation to the Hubbard model. 
The Fock term is absent in the Hubbard model. The solid line represents the undressed {U = 0) 
electron Green's function Qfj{iujn) and the dotted line represents the intrasite interaction U. The 
external legs just show the Matsubara frequency dependence and are not included in the analytic 
expressions. 
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FIG. 2. Magnetic phase boundaries for the second-order local approximation to the Hubbard 
model in three dimensions for two different values of the interaction strength U. The vertical axis 
shows the ratio of the temperature to the hopping constant. The horizontal axis shows the electron 
filling. The (dashed) solid lines denote the (in) commensurate transition from a paramagnetic to 
a spin-density- wave ground state in the thermodynamic limit. The (open) solid circles denote the 
corresponding transitions for a finite lattice of 80 unit cells for T > 0.05t and 500 unit cells for 
lower temperatures. The inset shows the corresponding z-component of the magnetic ordering 
wavevector q versus filling. The phase-diagram is symmetric around ng = 0.5. 
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FIG. 3. Magnetic phase boundaries for the first-order approximation to the Hubbard model in 
three dimensions for two different values of the interaction strength U. The vertical axis shows 
the ratio of the temperature to the hopping constant. The horizontal axis shows the electron 
filling. The (open) solid circles denote (in) commensurate transition from a paramagnetic to a 
spin-density-wave ground state for a finite lattice of 30 unit cells. The lines are a guide to the 
eye. The inset shows the corresponding z-component of the magnetic ordering wavevector q versus 
filling. The phase-diagram is symmetric around = 0.5. 
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FIG. 4. Comparison between the results of quantum Monte Carlo simulations^ (circles), the 
first-order (dashed line) and the second-order local approximation (solid line). Tcjt is the ratio 
of the transition temperature at half-filling to the hopping integral. U is the interaction strength. 
The second-order local approximation shows the same sign of the curvature as the QMC result, 
whereas the first-order curve shows the wrong sign of the curvature. Also, since the second-order 
approximation has the correct limiting behavior as ^ 0, this shows that the QMC data is 
overestimating Tc for small values of U. 
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FIG. 5. Comparison between the experimental data for chromium alloys and the results of 
the second-order local approximation. The (open) closed circles denote the experimental data for 
the (in)commensurate transitions in Cr alloys. The (dashed) solid lines denote the results of the 
second-order local approximation for the (in)commensurate transitions. The experimental data 
is best fit to the theoretical calculation when U = 5.5t. T/ is the transition temperature at the 
paramagnetic-commensurate-incommensurate phase boundary. Tj = S60K for experimental data 
and occurs when Cr is doped with 0.3 atomic percent of Mn. 
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